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The use of known analytic results for the continuum fermion determinants in QCD and 
QED as benchmarks for zero lattice spacing extrapolations of lattice fermion determi- 
nants is proposed. Specifically, they can be used as a check on the universality hypothesis 
relating the continuum limits of the naive, staggered and Wilson fermion determinants. 

1. Introduction 

Without internal quark loops such basic phenomena as color charge screening 
and fast quark fragmentation into hadrons cannot be formulated. Accordingly, the 
fermion determinant is one of the central quantities in any ab initio lattice QCD 
calculation that will confront experiment. 

It is proposed to use analytic results for fermion determinants in QCD and QED 
as a check on extrapolated lattice fermion determinants. In particular, the analytic 
results should be used to verify the universality hypothesis relating the continuum 
limits of the naive, staggered and Wilson fermion determinants for a single quark 
flavor: 1 

(det Wilson ) = det naive 

a— >0 

(detyviison) = det 

staggered j 

a— >0 

or stated somewhat differently, verify that 

detwilson) (det s taggered) i i (det nalve ) 16 = det ana lytic- 

a— >0 

It is understood that these are limits computed by taking ratios of the interacting 
and free determinants. Checks such as these would also be useful in establishing 
error bounds on determinant extrapolations to the continuum. One important point 
should be noted: det ana i y ti c is renormalization scheme dependant, and so it is nec- 
essary to have compatible lattice and continuum renormalization schemes when 
making the above comparisons. 

The analytic results, to be relevant and useful, should be free of ultraviolet and 
volume cutoffs. Since large quark mass can be dealt with by perturbation theory, 
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the results should be valid for all mass values or at least mass values small enough 
to connect with effective chiral theories. Finally, the background gauge fields should 
be as general as possible. Obviously these are severe constraints that greatly limit 
the available benchmark determinants. For lattice theorists such checks will require 
writing new code and convincing funding agencies that this is a worthwhile project. 

2. Defining the Fermion Determinant 

By fermion determinant we mean the ratio of determinants of the interacting and 
free Euclidean Dirac operators, det(f — gA + m)/dct(]f + m), defined by the renor- 
malized determinant on R d , namely 



det ren = exp(n 2 + n 3 + n 4 )det d +i(l - gSA), 



where 



lndet,j+i = Tr 



hi(l - 



E 

n=l 



(gSAy 



(i) 



(2) 



and S — + m) , 112,3.4 are the second-, third-, and fourth-order contributions 
to the one-loop effective action defined by some consistent regularization together 
with a coupling constant subtraction in II2 when d = 4; II4 is omitted when d < 4 
as well as II3 when d = 2. The operator SA is a bounded operator on the Hilbert 
space L 2 (R d , y/k 2 + m 2 d d k) for e r\ n>d L n (M. d ), in which case it belongs to the 
trace ideal C n for n > d [C n = {A|Tr(AU)"/2 < 00}]. This means that det^+i 
in (J2J can be expanded in terms of the discrete complex eigenvalues l/g n of the 
non-Hermitian compact operator SA: 



det d+1 (l-gSA)=H 



fjn 



exp 



d 

E 

fe=i 



(3) 



In QED, C-invariance and the reality of det rcn imply that the eigenvalues of SA 
appear in quartets ±l/g n , ±l/gr*; in addition, II3 is absent in qTp . It is always 
assumed that m ^ 0. The m = case raises nontrivial problems, and most of the 
above statements, including are false in this case. 2 The reader is referred to 
Refs. 2-4 for a review of these basic results. 

Definition Q of det ron gives back what one gets using the Feynman rules, and 
so det r en is what physicists mean by a determinant. Closely related to Q is the 
Schwinger proper time representation of the determinant provided it has a coupling 
constant renormalization subtraction included. These are to be contrasted with the 
zeta function definition whose field-theoretic meaning is not evident. The modified 
Fredholm determinant m the definition of det ren is necessary because the loop 
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diagrams representing Tt(SA) 71 for n = l,..,d are not denned by themselves. We 
have summarized this basic material because some lattice determinants in the lit- 
erature, such as det(f — gA + m) = [] n (A n + m), make no sense when extrapolated 
to zero lattice spacing. 

3. Analytic Results 

Because det ren has no zeros on the real g-axis and det rcn (g = 0) = 1, it is positive 
for real g. By inspection it is an entire function of <?, and because SA € C n , n > d, 
it is of order d. This means that for suitable constants A(e), K(e) and any complex 
value of g, |det rcn | < A(e) exp (K(e)\g\ d+,i ) for any e > 0. For d = 2,3 both the 
QCD and QED determinants satisfy det rcn < 1 for real g and all values of m > 0. a 
This is known as the " diamagnetic" bound 3-6 which can be interpreted in QED as a 
reflection of the paramagnetism of charged fermions in a magnetic field. Therefore, 
in these cases det ren definitely does not achieve its maximal growth on the real 
<?-axis. Already at this point we have exhausted all that can generally be said so 
far about fermion determinants in QCD and QED. 

Keeping to the criteria we have stated at the end of Sec. ^ we have only two 
cases to report. The first deals with QCD's fermion determinant in the instanton 
background: 



The t] symbols arc defined in the Appendix of Ref. 7; p is a scale parameter. For 
mp << 1, 

lndet rcn = ln(mp) — - ln(/ip) + 2a (~ \ + [\n(mp) + 7 — In 2] (mp) 2 + 0(mp)' 1 , (5) 
where 



3 ^" \2 



i)=i( 7 + m*)-M_g = 0.145873, (6) 

7 is Euler's constant and /x is the subtraction point of the coupling constant g. The 
first three terms in JSJ were obtained by 't Hooft. 7 Later, Carlitz and Creamer 8 
obtained the (mp) 2 ln(mp) term. Then Kwon, Lee and Min 9 verified the calculation 
in Ref. 8 and also obtained the remaining terms of order (mp) 2 . Note that the 
coefficient of the leading ln(mp) term, 1, is the chiral anomaly in this case, which 



a This bound is initially derived on a lattice. There remains the technical problem of proving that 
the zero lattice spacing limit of the lattice fermion determinant coincides with det rcn defined by 
til and This is partially dealt with in Refs 3,4. A continuum proof of detren < 1 with g real 
for d = 2, 3 is needed. 
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should be of special interest in extrapolations of lattice calculations of the instanton 
determinant. 

Recently, Dunne, Hur, Lee and Min 10 have proposed a partial-wave WKB phase- 
shift method for calculating the instanton determinant for all values of m. Although 
their predicted value for a(l/2) is 0.137827 to third order in the WKB approxima- 
tion, inclusion of higher order terms may improve this result. If so this would raise 
the question of why the WKB analysis should work for other than large values of 
m. 

The inclusion of the subscript ren on the determinant in Q is provisional 
because it was not calculated following (JJJ and (J2J. Instead, the 2a(l/2) term in 
lO was obtained using Pauli-Villars regulator fields b with space-time dependent 
masses followed by a change in the subtraction point to the constant mass /i in JSJ . 
It would be reassuring if the instanton determinant were recalculated for m =/= 
using Q and (J2J) or the Schwinger proper time representation without space-time 
dependant masses and verify that the 2a(l/2) term emerges when rap << 1. 

It is easy to generalise the first two terms in (0 to any square- integrable self-dual 
background gauge field with range p interacting with a massive fermion. Referring 
to Eqs. (3.5) and (3.10) in Ref. 12 one sees by inspection that for QCD 4 and QED 4 



J^o m2 & lndet ™ = 32^ - JLj j (7) 

The last term in J7J) arises when renormalization is carried out by a momentum 
space subtraction at p 2 — 0. Subtracting instead at p 2 — p 2 gives for mp « 1 



\ndet rcn = — [ln(mp) - - hx(pp)} J d*x trF 2 „ + i?(m 2 , p 2 ), (8) 

where R(0,p 2 ) is a p, and p— independant constant that depends of F^ v . 

The second analytic result available is the fermion determinant in QED 2 for the 
background field = e^d^ip, where ei2 = —£21 = 1 and B = —d 2 (p. The magnetic 
field B is arbitrary except that B(r) > with continuous first and second derivatives 
and B(r) — for r > a with J Q ° dr rB 2 (r) < 00. In polar coordinates, A r = 0, 
A g = $(r)/27rr, where $(r) is the total flux of B out to r. Hence Ag = $/2irr, 
r > a, where $ is £Ts total flux. Then for the sequence of limits indicated we obtain 



b Diffcrcnt regulator schemes can affect the 2a(l/2) term in J^J. If dimensional renormalization is 
used instead then 

-(la4ir-7) = -0.35952. 
Eqs (13.7), (13.9) and (13.12) in Ref. 7 relating to this change are incorrect. 11 









-G) 
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the following result: 13 

Urn lim lndet rcn = -1^1 In ( -J^L ) + 0(|e4>|, (ma) 2 |e$| ln(|e$|)). (9) 

|e$|>>lma<<l An \(ma) 2 J 

The minus sign is in accord with the diamagnetic bound cited in Sec.[5J and the 
coefficient of the ln(ma) term, e$/27r, is the two-dimensional chiral anomaly, assum- 
ing e > 0. The analysis in Ref. 13 is not sufficient to rule out the (ma) 2 1 e<E» | ln|e$| 
term in the remainder in 0; it may only be (ma) 2 |e$|. 

This result is therefore valid for a whole class of background gauge fields in the 
nonperturbative region of small mass and strong coupling. Note that the II2 term 
in is cancelled by det% in Q and (J2J). We believe this is a general result and 
that II4 will be cancelled by det$ in d = 4. 

Finally, although result 10 is for an Abelian background field in d = 2 it remains 
a challenge for lattice QCD fermion determinant algorithms - which we take to 
include the discretization of the Dirac operator adopted - to reproduce it. 

Added Note: The authors cited in Ref. 10 have substantially improved their 
result in hep-th/0410190 ("Precise Quark Mass Dependance of the Instanton De- 
terminant"). 
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